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COMMON-PATH HETERODYNE LASER-INDUCED THERMAL ACOUSTICS 
FOR SEEDLESS LASER VELOCIMETRY 


ROGER C. HART*, G. C. HERRING*, and R. JEFFREY BALLA* 

Abstract. We demonstrate the use of a novel technique for the detection of heterodyne laser-induced 
thermal acoustics signals, which allows the construction of a highly stable seedless laser velocimeter. A common- 
path configuration is combined with quadrature detection to provide flow direction, greatly improve robustness to 
misalignment and vibration, and give reliable velocity measurement at low flow velocities. Comparison with Pitot 
tube measurements in the freestream of a windtunnel shows root-mean-square errors of 0.67 m/s over the velocity 
range 0-55 m/s. 

Key words. LIT A, laser velocimetry, laser anemometry, laser-induced thermal acoustics 

Subject classification. Experimental Fluid Dynamics 

L Introduction. Non-intrusive optical methods for measuring fluid flow velocity are of great importance 
to the experimental fluid dynamics and aerodynamics communities. However, the only optical velocimetry methods 
to have found wide use, laser doppler velocimetry [ 1 ] and particle image velocimetry [2], require the introduction of 
small (~ I pm diameter) seed particles into the flow to serve as light scatterers. Seeding is not feasible in some wind 
tunnels due to concerns over removal of spent seed, clogging of flow-straightening screens, or abrasion of finely 
polished surfaces. Additionally, there are regions in airflows of interest, such as vortices over delta wings or behind 
leading-edge slats, where useful seed concentrations are difficult or impossible to achieve. Thus a need exists for a 
seedless laser velocimetry method. Previous work has demonstrated the potential of laser-induced thermal acoustics 
(LITA) to fill this need [3-6]. However, none of the approaches employed so far for heterodyne detection appears 
suitable for routine use in environmentally adverse windtunnel environments by workers who are not highly skilled 
optics researchers. Here we demonstrate a novel means of implementing LITA velocimetry that has the intrinsic 
stability and robustness to allow ihe construction of a useful, fieldable instrument. 

2. Theory of LITA Signal Detection. LITA is a pump-probe process: a pump laser creates a periodic 
perturbation in the medium, which serves to diffract the beam of a probe laser to a detector. In our apparatus the 
beam of a Q-switched Nd:YAG laser operating at wavelength X = 1064 nm is split by a 50/50 beamsplitter to 

produce two pump beams which are made to focus and cross at a single point at angle 20 = 1.4 by a lens (Fig. 2.1). 
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Interference fringes of period A = X /? „, w/ ,/2sin0 are formed at the beam intersection. As the laser wavelength is not 
resonant with any absorptive transition in the medium (air) the only optical interaction is electrostriction, which 
creates two counter-propagating acoustic waves or gratings traveling in the ±x direction at speed of sound v s with 
respect to the medium. The acoustic wavelength is equal to the optical fringe spacing and the frequency 
f B = v s / A is about 7.8 MHz at room temperature. The beam of the CW probe laser at wavelength X /wAt , = 532 nm 



Fk.. 2.1. Beam crossing geometry and schematic of detection scheme. A: attenuator; B: beamstop; C: Focal/detector plane; 
LO: local oscillator beam. 


is incident on the waves at the Bragg or phase-matching angle <|> where A = X he j2sm^ . The signal beam is 
diffracted at angle <j) (Fig. 2. 1 ) with Doppler shifts Aco^ = — q • (v tw| - v s ) and Aco_ = -q • (v ton + v 5 ) due to the 
motion of the gratings, w'here v con is the flow velocity, v s = v s x , and q = qx is the acoustic grating vector with 
q -2n/ A . 

2.1. Collinear Local Oscillator Method. The simplest implementation of LITA velocimetry (not shown) 
requires the introduction of a local oscillator (LO) beam collinear with the signal with the same frequency to as the 
probe. The field is 

(2.1) A{t) = ^ /0 exp(/'co/)-f A s e p, (exp(/(a> + Ao>_ /)) + exp(/(to + Ao ) 4 /)))+c.c. , 





where A lo and A s are the real amplitudes respectively of the LO and signal beams, P = a v s with a the acoustic 
amplitude absorption coefficient [7], and c.c. means complex conjugate. Averaging over a time long compared to 
the optical period but short compared to 1 if B and integrating over the area of the beams, the lime dependence of the 
detected power is 


(2.2) P(t) = Aj () + 2A U) A s e ^ w [cos(Aco_ /) + cos(Aco_ /)] , 

where we assume A lo » A s for shot-noise-limited detection and so neglect the term x A$ . Neglecting also the 
DC term, the real part P R (co) of the Fourier transform of P(t) is a sum of four Lorentzians 


p _ + _ p 

P“+(eo±Aco )" p “ + (d) ± A(o.. )“ 

with peaks at ±Aco^ and±Aco . Although the Doppler shifts are signed quantities, the detected heterodyne signal 
contains information only on the absolute values. Analysis (e.g. discrete Fourier transform or time domain fitting 
via Levenburg-Marquardt or Prony’s method [6]) of real-valued data described by Eq. 2.2 thus allows independent 
determination only of |Aco_ | and |Aco . | , yielding v s and |vy | , where v, = q * v aw is the component of flow velocity 
along q . 

Laboratory evaluation of this collinear LO method of implementing LITA velocimetry [6] revealed the 
following disadvantages: 1 ) proper alignment of the LO and signal beams is difficult to achieve and maintain; 2 ) the 
apparatus is very sensitive to vibration since the signal and LO beam must follow different paths on different optics; 
3) no direct means of determining flow' direction is available since only the absolute values of Aw and Ato are 
available as shown above; 4) both accuracy and precision are very poor for flow' velocities less than -20 m/s due to 
the difficulty of accurately discriminating the nearly degenerate frequencies |Aco | and | Aco ^ | in the presence of 

noise. Although various technical means of ameliorating these difficulties are available, this approach seemed 
unsuitable to be the basis of a robust instrument with a useful velocity dynamic range. 

2.2. Grating Demodulation with Quadrature Detection. A solution to these four problems was found 
with two modifications: a non-collinear geometry, and quadrature detection. A lens recrosses and refocuses the 
signal and attenuated probe beams in plane C (Fig. 2. L). Assuming for simplicity unit magnification, the field at C 
is 

(2.4) A ( /) = A l() exp ( / (k U) * r - o> t ))+ 

A s e ^ ’ exp (/' (yf v • r - (co + A co V ))+ 

A s e ^ 1 exp (/ ( A' v • r - (co + A co , )f))+ c.c , 


(2.3) 


P R (<5) = 2A LO A s 
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where k LO zvi<\ k s are the wavevectors for LO and signal beams with = 2n/k prohe . With | k LO - = q (due 

to the phase matching between pump and probe), and neglecting the term <x a\ , the intensity in plane C is then 

(2.5) /(*,/) = A~ U) + 2 A lo A$e ^ 1 [cos(^v + Aco /)+ cos(^.v + Au>^/)] , 

which describes interference fringes of spatial period A traveling at speeds Aco + /r/ = v F - v v and 
Aco /q = v F + v v ; that is, the lens produces an image of the acoustic waves. No modulation at frequencies 
Ato_ and Aco^ would be seen by a broad-area detector (width » A ) placed at C. However, if a series of slits (e.g. 

a Ronchi ruling) parallel to the fringes with period A -2D is placed on the surface of the detector (Fig. 2.2) the 
modulation is recovered as bright fringes pass alternately across transparent and opaque regions. We refer to this 
detection scheme as grating demodulation. Taking advantage of the equal periodicities of the fringes and the 
Ronchi ruling, the time dependence of the detected power P(t ) may be written 

(2.6) P(t) = F r°Y (x,t)dx 

where F is a constant accounting for the number of fringes and the height of the beams, and a specifies the 
transverse location (phase) of the ruling with respect to an arbitrary origin (Fig. 2.2). After some manipulation one 
finds 

(2.7) P(t) = hAl 0 

- 2cA LO A s e~^ t (sin((j))cos(Aa)_/)+ cos(cp)sin(Aco_/)) 

-2 cA l0 A s e~^ * (sin(q) )cos(Aco + /) + cos(tp )sin(Aco + /)), 
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where cp = na/ D and b and c are constants. Thus by shifting the ruling transverse to the optical axis (changing a) 
the phase of the modulation in the signal may be varied; this is equivalent to shifting the phase of the heal 
oscillator . This provides a convenient means for implementing dual-channel phase-sensitive detection with two 

phases of the LO differing by 90°, known as quadrature detection. If the converging signal and LO beams (Fig. 1) 
are divided by a suitable beamsplitter so that two separate beam intersections are formed on two separate rulings and 
detectors, then by setting the relative transverse shift between the rulings to A/4 both quadrature components of the 

signal may be recovered simultaneously. That is, for channel />,(/) let a = 3D/2 in Eq. 2 and for channel i %( t) 
a = D , giving 

(2-8) P {] (t) = bAl a + 2cA LO A s e p '(cos(Aco /) + cos(Ato ./)), 

Pi)o(t) = bAj 0 + 2cA LO A s e~^' (sin(Aco f) + sin(Aco^)). 

Forming the complex signal Z(t) = P f) (t) + iP 90 (t ) , we find the real part of the Fourier transform now consists of 
only two Lorentzians 


(2.9) 


Z^oT) = 2 cA L() A s 


P 2 + (co - Aco ) 2 p 2 + (co - Aco^ ) 2 


with peaks at Aoj and Aco^ . Thus analysis of Z(t ) can recover the sign as well as magnitude of the Doppler shifts 
(Fig. 2.3). For v F = 0 the spectral peaks occur at (S = ±qv s =±w B ^ ±2nf B ; as v F varies the peaks w ill shift 
simultaneously towards more positive or more negative frequency (depending on the flow direction) but will 
maintain the separation 2co^ . Flow' and sound velocities are found as v F =(Aco 4 +Ato )/2 q and 

v s = ( Aco » " A(0 )/2<? • Note that Mach number M = v F /v s is independent of q , so if flow temperature and thus 
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v 5 are known v F may be found without calibration of A . 

This grating demodulation/quadrature detection scheme remedies the defects found in the collinear LO 
method. Since signal and LO beams are incident on the same optics (common path), detection is intrinsically stable 
against both misalignment and vibration. Flow direction is found directly. Resolution at low velocities is greatly 
improved since the frequencies Aoo and Aco t are no longer nearly degenerate. 

3. Experimental Results. A compact velocimeter employing this detection scheme was used at NASA 
Langley Research Center’s Basic Aerodynamics Research Tunnel to map the flow behind a rearward-facing step [8]. 
As part of this program, multiple comparisons of freestream velocity V L i T a (average of -350 shots) measured by 
LITA with freestream velocity V Pi(ot measured by Pitot tube were acquired (Fig. 3.1). Defining error as V Pjtor V LITA , 
the root-mean-square error (RMSE) of this data set is 0.67 m/s. For V P j, 01 > 30 m/s and for V P , 101 = 0, the RMSE is 
< 0.4 m/s, while for intermediate velocities RMSE is 0.84 m/s. Note that a velocity of 1 m/s corresponds here to a 
Doppler shift of 23 kHz in a carrier frequency of 5.6 x 10 14 Hz. 

We have demonstrated an elegant, intrinsically stable means of performing heterodyne LITA velocimetry with 
quadrature detection that will allow the construction of a robust seedless laser velocimeter. 
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Fi(i. 3.1. Error Vp, u ,r Vuta for measurements in windturmel freestream. Error bars show one standard deviation of the mean. 
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